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摘要 :【 目 的】 本 研究 拟 选 择 适 合用 于 分 析 松 墨 天 牛 Monochamus alternatus 化 学 感受 组 织 中 基因 表 
达 的 内 参 基因 。[【 方 法 ] 依 据 转 录 组 测序 结果 进行 内 参 基因 鉴定 ,利用 RT-qPCR 技术 分 析 内 参 基 因 
在 松 墨 天 牛 不 同 发 育 阶 段 和 不 同性 别 化 学 感受 组 织 间 的 表达 差异 ,并 利用 软件 geNorm, 
NormFinder 和 BestKeeper 比较 其 表达 的 稳定 性 。 【结果 】 松 墨 天 牛 转录 组 中 鉴定 出 9 个 候选 内 参 
基因 (Actin, TUB, 18S rRNA, RPS27A, RPS3, RPL10, AK, GAPDH 和 EFIA) ,其 中 后 7 个 候选 内 
参 基因 在 松 墨 天 牛 中 被 首次 鉴定 , 松 墨 天 牛 候选 内 参 基因 和 其 他 昆虫 相应 基因 的 同 源 性 很 高 。9 
个 候选 内 参 基因 引物 均 具有 良好 的 扩 增 效率 ,18S rRNA. 的 表达 水 平 最 高 ,EF14 的 表达 水 平 最 低 ; 
18S rRNA £e Actin 在 不 同样 品 间 的 表达 水 平 差异 最 大 ,C4PD fa TUB 表达 水 平 在 不 同样 品 间 差 异 
Jh. geNorm 和 NormFinder 软件 分 析 认 为 ,GAPDH 是 最 稳定 的 内 参 基因 ,TUB 是 较为 稳定 的 内 参 
基因 ,18S rRNA 和 Actin 是 最 不 稳定 的 内 参 基因 ;BestKeeper 软件 分 析 认 为 ,GA4PDH 和 TUB 是 合适 
的 内 参 基 因 ,18S rRNA 和 Actin 是 不 适合 的 内 参 基因 。 最 适合 校正 松 黑 天 牛 化 学 感受 组 织 中 基因 
表达 数据 的 内 参 基 因数 量 为 2 个 , 即 G4PDH 和 TUB, 并 且 这 样 的 内 参 基因 组 合 可 以 用 于 不 同 发 育 
阶段 和 不 同性 别 的 不 同化 学 感受 组 织 。[【 结论 】] 本 研究 结果 为 利用 RT-qPCR 技术 准确 分 析 松 墨 天 
牛 和 其 他 天 牛 基 因 包 括 化 学 感受 组 织 基因 相对 表达 量 的 内 参 基因 选择 提供 了 重要 参考 。 
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Identification and selection of valid reference genes for assaying gene 





expression in the chemosensory tissues of Monochamus alternatus 


( Coleoptera: Cerambycidae) by RT-qPCR 

FENG Bo, GUO Qian-Shuang, MAO Bi-Peng, DU Yong-Jun ( Institute of Health and Environmental 
Ecology, Wenzhou Medical University, Wenzhou, Zhejiang 325035, China) 

Abstract: [ Aim] This study aims to select suitable internal control genes for gene expression analysis in 
chemosensory tissues of Monochamus alternatus. [ Methods] Candidate reference genes were identified 
from our transcriptome results, and their expression patterns in chemosensory tissues were investigated by 
RT-qPCR in different developmental stages and both genders. Their expression stabilities were also 
compared with each other by using softwares geNorm, NormFinder and BestKeeper. [Results] Nine 
candidate reference genes were identified in M. alternatus, including Actin, TUB, 18S rRNA, RPS27A, 
RPS3, RPL10, AK and GAPDH, of which the last seven were firstly identified in this insect species. 
Candidate reference genes in M. alternatus have high nucleotide sequence identity with those from other 


insects. Designed primers of nine candidate reference genes worked well with high amplification efficiency 
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and specificity to M. alternatus. It was found that 18S rRNA had the highest expression level while EFIA 
had the lowest expression level in M. alternatus. 188 rRNA and Actin had the most variable expression 
levels between studied chemosensory tissues, but GAPDH and TUB had the least variable expression 
levels. GAPDH was the most stable reference gene, TUB was a moderately stable reference gene, and 
18S rRNA and Actin were the least stable reference genes in different chemosensory tissues, which was 
referred from analysis with softwares geNorm and NormFinder. Analysis using BestKeeper software showed 
that GAPDH and. TUB. but not 18S rRNA and Actin were the suitable reference genes in chemosensory 
tissues of M. alternatus. The combination of two reference genes, GAPDH and TUB , was referred as the 
best internal control combination genes for normalization of RT-qPCR data conducted with chemosensory 
tissues of different developmental stages (larva, pupa and adult) and different genders of M. alternatus. 
[Conclusion] The results of this study provide important data for selecting internal reference genes for 


gene expression analysis of M. alternatus and other longhorn beetles, and also for gene expression 


analysis of chemosensory tissues of other insects by RT-qPCR. 


Key words: Monochamus alternatus; reference gene; chemosensation; development; RT-qPCR 


实时 荧光 定量 反 转 录 PCR(real-time quantitative 
reverse transcription-PCR ，RT-qPCR ) 是 分 子 生物 学 
领域 研究 基因 表达 的 最 有 效 方法 之 一 (Yan et al., 
2012) 。 在 RT-qPCR 实验 过 程 中 ,起 始 RNA 质量 、 
cDNA 合成 效率 等 方面 的 差异 可 能 影响 结果 ,为 了 
去 除 这 些 影 响 因 素 ,需要 引入 内 参 基 因 对 数据 进行 
校正 和 标准 化 处 理 ( Chen et al., 2012) 。 选 择 适 当 
的 内 参 基因 对 于 分 析 结 果 的 准确 性 至 关 重 要 ,不 合 
适 的 内 参 基 因 可 能 导致 错误 结论 (Jian et al., 
2008) 。 因 此 在 选择 内 参 基 因 之 前 ,需要 对 候选 内 
参 基 因 的 有 效 性 进行 评 佑 。 

理想 的 内 参 基 因应 该 稳定 表达 于 不 同类 型 的 组 
织 和 同一 组 织 的 不 同 处 理 中 ,其 表达 水 平 不 受 任何 
内 源 性 或 外 源 性 因素 的 影响 (Janska et al., 2013) 。 
过 去 普遍 认为 ,维持 细胞 正常 生命 代谢 所 必需 的 持 
家 基因 在 不 同 组 织 及 不 同 实验 条 件 下 都 具有 稳定 的 
表达 模式 ,是 理想 的 内 参 基因 (Thellin et al., 1999) 。 
但 越 来 越 多 的 研究 表明 ,这 些 内 参 基因 在 不 同 物种 、 
不 同 发 育 阶 段 .不 同 组 织 及 不 同 试验 条 件 下 的 表达 
稳定 性 存在 差异 (de Jonge et al., 2007; Sun et al., 
2009; Huis et al., 2010) 。 研 究 者 们 开始 筛选 特定 
物种 体内 稳定 表达 的 内 参 基 因 。 

已 经 报道 过 内 参 基 因 筛 选 的 物种 包括 小 沫 邮 
Plutella xylostella (Fu et al., 2013), 棉铃 IE 
Helicoverpa armigera ( Shakeel et al., 2015) ,甜菜 夜 
IR Spodoptera exigua (Zhu et al., 2014) 和 其 他 蛾 类 
(Teng et al., 2012) , E X JR =Œ IH FH. Diabrotica 
virgifera ( Rodrigues et al., 2014), K Js np FH 
Colaphellus bowringi ( Tan et al., 2015) , 5 4$ 3E FH IK 
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Leptinotarsa decemlineata (Shi et al., 2013) 等 甲虫 及 
3 UNIES Mylabris cichorii ( Wang et al., 2014) 等 ， 
到 目前 为 止 还 没有 关于 天 牛 内 参 基 因 筛 选 的 报道 。 
在 进行 内 参 基 因 筛 选 时 ,最 党 使 用 的 组 织 为 消化 道 
(Fu et al., 2012; Peng et al., 2012; Lu et al., 2013) 
和 整个 虫 体 (Lord et al., 2010; Rodrigues et al., 
2014; Chen and Lu, 2014) ,也 包括 头 部 (Teng et al., 
2012) 脂肪 体 (Peng et al., 2012; Shi et al., 2013) , 
Ub EE (Peng et al., 2012) ftl S£21 细胞 (Salem et al., 
2014) , ELLAS B 5 ( Teng et al., 2012) Ji && ( Teng et 
al., 2012) 表皮 (Shi et al., 2013) 、 丝 腺 (Teng et al., 
2012) 、 血 淋巴 (Lu et al., 2013; Zhu et al., 2014) 等 
组 织 , 但 是 到 目前 为 止 几乎 没有 研究 涉及 到 化 学 感 
受 组 织 。 

H&K Monochamus alternatus Hope 是 国际 上 
公认 的 重要 林业 有 害 生 物 松 材 线虫 Bursaphelenchus 
xylophilus 的 主要 传播 媒介 ,也 是 一 种 严重 危害 林木 
的 性 干 害 虫 (Hu et al., 2013) 。 串 觉 在 松 墨 天 牛 的 
取 食 、 产 卵 、 寄 主 转移 等 个 体 行 为 中 具有 重要 的 行 》 
调控 作用 ,寄主 植物 释放 的 蓝 烦 类 挥发 性 次 生物 质 
能 够 引诱 松 黑 天 牛 成 虫 取 食 (Ikeda et al., 1980) , 生 
理 衰弱 和 新 死 松树 释放 出 来 的 昔 燃 化合 物 和 乙醇 引 
VETASECK ^F WEE rh y^ B ( Hao et al., 2006) , 松 黑 天 牛 
成 虫 可 能 通过 雄性 个 体 分 泌 的 聚集 信息 素 进行 通信 
识别 (Teale et al., 2011)。 本 人 研究 利用 RT-qPCR 方 
法 检测 松 墨 天 牛 不 同 发 育 阶段 和 不 同性 别 的 化 学 感 
受 组 织 中 候选 内 参 基因 的 表达 水 平 ,利用 生物 学 软 
TF geNorm, NormFinder 和 BestKeeper 评估 候选 内 参 
基因 的 稳定 性 ,筛选 出 适合 于 松 墨 天 牛 化 学 感受 组 
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织 中 定量 表达 分 析 的 内 参 基因 。 cDNA 于 -20 % 保存 ,用 于 RT-qPCR。 


l 材料 和 方法 


1.1 供 试 昆虫 

松 墨 天 牛 师 采 自 浙 江宁 波 松树 ,将 幼虫 臂 出 之 
后 放 入 5 mL 指 形 管 中 ,在 温州 医科 大 学 健康 与 环境 
生态 研究 所 人 工 气候 室 进行 饲养 。 幼 虫 喂 以 松树 碎 
届 。 饲 养 温度 为 26 +10 ,相对 湿度 为 75% x596, 
光 周 期 为 14L: 10D。 成 虫 羽化 后 单 尖 饲养 在 50 mL 
指 形 管 中 , 喂 以 新 鲜 松 校 ,每 3 d 更 换 一 次 松 枝 。 收 
READE ALMS d 的 肾 及 取 食 10 d 的 成 虫 用 于 
实验 。 

在 显微镜 下 ,利用 显 微 解剖 和 剪刀 分 别 取 下 幼虫 、 
AA E RR A ARE R R KO fik fA A T SE CE E , ME JL R 
和 雄 成 虫 的 足 , 雄 成 虫 的 头 部 和 腹部 末端 等 化 学 感 
受 组 织 。 在 液 氮 中 快速 冷冻 之 后 ,样品 放 入 - 80€ 
冰箱 保存 。 
1.2 总 RNA 的 提取 和 cDNA 合成 

采用 试剂 TRNzol-A * (天 根 生化 科技 有 限 公司 ， 
中 国 ) 进 行 总 RNA 的 提取 ,具体 方法 如 下 :在 液 毛 
中 将 10 mg 组 织 研 碎 , 加 入 1 mL TRNzol-A -试剂 , 静 
5 min 后 ,加 入 0.2 mL 和 氯仿 ,4% 12 000 r/min 离 
心 15 min 后 ,向 上 清 液 中 再 次 加 入 0.2 mL 氯仿 ， 
4*C 12 000 r/min 离心 15 min 后 ,向 上 清 液 中 加 入 
等 体积 异 两 醇 , -20% 放置 2h 后 ,4% 12 000 r/min 
离心 15 min ,利用 1 mL 7596 乙醇 洗涤 沉淀 2 次 ,最 
后 加 入 30 ~ 100 uL 去 RNase 的 蒸馏 水 溶解 RNA。 
之 后 用 1% 琼脂 糖 凝 胶 电 泳 ,检测 总 RNA 完整 性 。 
并 用 紫外 分 光 光 度 计 NanoDrop2000 (ThermoFisher， 
America) 进行 浓度 和 纯度 分 析 。 为 了 保证 实验 的 准 
确 性 ,电泳 结果 中 均 可 见 清晰 的 2 条 带 (28S rRNA 
和 18S rRNA) ,其 条 傍 亮 度 比值 均 大 于 3: 1 ,紫外 分 
JOGE OD ODs f fE 1. 70 ~2.05 之 间 ， 
0D:oZOD2o 值 大 于 1 的 样品 用 于 后 续 实 验 。 用 
DNase | (NEB, America) 处 理 去 除 RNA 中 的 基因 
组 DNA 污染 。 

利用 逆转 录 试 剂 盒 ReverTra Ace qPCR RT Kit 
(Toyobo, Japan) 按照 使 用 说 明 书 进行 RNA 的 逆转 
录 , 上 有 具体 方法 如 下 :将 0.5 pg 总 RNA 在 65 伟 条件 下 
处 理 5 min 后 , 置 于 冰 上 ,分 别 加 入 5 x RT 缓冲 液 2 
HL, 逆转 录 混 合 物 0.5 uL, 引物 混 合 物 0.5 uL, Bf 
利用 去 RNA 酶 的 水 补足 10 pL, 在 37 今 条件 下 逆转 
录 15 min, 最 后 在 98 人 条 件 下 处 理 5 min 使 酶 失 活 。 


































































































1.3 引物 设计 和 合成 

从 我 们 的 转录 组 测序 结果 中 ,选择 了 9 个 常用 
的 内 参 基 因 :a- 微 管 蛋白 基因 (alpha-tubulin gene, 
TUBA) , i2 & -ECWEUK SR H S27a 基因 (ubiquitin-40S 
ribosomal protein S27a gene, RPS27A4) ,60S 核糖 体 
蛋白 Ll0a 基因 (60S ribosomal protein L10a gene, 
RPL10) , 精 氮 酸 激酶 基因 (arginine kinase gene, 
AK) ,甘油 醛 磷 酸 脱 氧 酶 基因 (glyceraldehyde-3- 
phosphate gene, GAPDH), 延伸 子 1A X 
C elongation factor 1 alpha gene, EFIA) ,40S 核糖 体 
蛋白 SS X (40S ribosomal protein S3 gene, 
RPS3 ) ,18S 核糖 体 RNA (18S rRNA ) 和 肌 动 蛋白 基 
Pd (Actin) 。 根 据 得 到 的 基因 序列 ,利用 软件 Primer 
Premier 6.0 ( Premier, Canada) 设计 定量 分 析 引 物 
( 表 2) ,各 引物 的 设计 长 度 为 20 ~ 22 个 碱 基 , 退 火 
温度 为 60%C , 扩 增 产物 长 度 小 于 250 pp。 引 物 由 上 
海 桑 尼 生物 技术 有 限 公司 合 成 。 
1.4 RT-qPCR 

使 用 TOYOBO 公司 的 Thunderbird SYBR qPCR 
Mix 试剂 盒 进行 RT-qPCR。PCR 反应 体系 包括 :2 x 
Thunderbird SYBR qPCR Mix 10 pL,4 4% f TE HJ 
cDNA 模板 1 pL; 上 下 游 引物 (10 pmol/L) 各 1pL， 
无 核酸 酶 的 水 补 齐 到 20 uL, PCR 反应 在 Bio-Rad 
公司 的 C1000"" Thermal Cycler PCR 仪 上 完成 。 
PCR 反应 程序 为 :95% 预 变性 3 min;95Y 变性 10 s, 
60Y 退火 和 延伸 共 30 s, 在 延伸 阶段 检测 欧 光 强度 ， 
收集 信号 ,40 个 循环 。 扩 增 结束 后 ,做 熔 解 曲线 分 
析 以 检测 扩 增 产物 的 特异 性 ,温度 为 60 ~ 95C , 升 
温 速 度 为 0.5%CAS。 熔 解 曲线 为 单 峰 时 用 于 后 续 分 
析 。 每 次 实验 设置 3 个 平行 重复 孔 。 幼 虫 . 师 、 雌 成 
虫 和 雄 成虫 的 触角 和 下 展 ( 下 上 5) y ME R AAE R E 
的 足 , 雄 成 虫 的 头 部 和 腹部 末端 共 12 个 样品 ,每 个 
样品 进行 3 次 生物 学 重复 。 
1.5 标准 曲线 的 绘制 和 引物 扩 增 效率 验证 

将 反 转 录 得 到 的 cDNA 模板 稀释 2 倍 后 ,依次 
稀释 5 个 梯度 ,每 个 梯度 稀释 5 倍 , 分 别 得 到 5 CI 
反 转 录 产 物 的 2 倍 稀释 液 ), 57, 5”, 5“ 和 5“ 共 
5 个 浓度 的 cDNA 模板 (Su et al., 2013) ,按照 1.4 
节 进 行 RT-qPCR。 依 据 结果 ,利用 线性 相关 方法 分 
析 cDNA 模板 浓度 和 Ct 值 之 间 的 线性 关系 ,得 到 相 
关系 数 ( 尼 ) 和 和 斜率 (slope), 利用 公式 = 
(10 7^" -1) x100 计算 引物 的 扩 增 效率 (EE) 
( Radonic et al., 2004) 。 
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1.6 数据 分 析 

根据 标准 曲线 将 RT-qPCR 得 到 的 Ct 值 转化 为 
相对 定量 数据 ,利用 geNorm ( Vandesompele et al., 
2002) ,BestKeeper ( Pfaffl et al., 2004) 和 NormFinder 
( Andersen et al., 2004) 软件 严格 按照 软件 使 用 说 
明 对 各 候选 内 参 基因 进行 筛选 。 在 利用 geNorm 分 
析 时 ,表达 稳定 值 M 小 于 0.5 时 即 表示 该 基因 为 稳 
定 表 达 ( Hellemans et al., 2007) ) ,配对 差异 值 小 于 
0.15 时 表示 不 需要 再 增加 内 参 基 因 的 数量 
(Vandesompele et al., 2002) 。 在 利用 BestKeeper 软 
件 分 析 时 , 当 标 准 变异 系数 (5D) 大 于 1, 或 者 P 值 
大 于 0.05 时 , 则 该 基因 不 适合 作为 内 参 基因 (Pfaff 
et al., 2004) 。 不 同 发 育 阶段 分 析 时 选用 松 墨 天 牛 幼 
虫 ` 雄 师 和 雄 成 虫 的 触角 和 下 层 须 的 扩 增 数据 ,不 同 
性 别 分 析 时 选用 上 肉 雄 成 虫 触角 、 下 层 须 和 足 的 扩 增 数 
据 , 不 同 组 织 分 析 时 选用 雄 成 虫 触角 、 足 、 下 层 须 、 腹 
部 末端 和 头 部 的 扩 增 数据 ,筛选 出 最 优 的 内 参 基 因 。 









































2 结果 


2.1 BÆRES 个 内 参 基因 分 析 

根据 我 们 前 期 转录 组 的 测序 结果 , 共 筛 选 到 9 
个 内 参 基 因 ( 表 1) , 除 TUB 和 Actin 外 ,其 余 7 个 基 
因 在 松 墨 天 牛 中 首次 被 鉴定 ,各 基因 在 CenBank 中 
的 登录 号 见 表 1o RPS27A, RPLIO, AK, GAPDH, 
EFI1A 和 RPS3 的 基因 长 度 分 别 为 632, 754, 1 475, 
1 335, 2 897 和 832 bp ,包含 长 度 分 别 为 471，654 ， 
1 071, 1 002, 1 773 和 729 bp 的 完整 开放 阅读 框 ， 
分 别 编码 156, 217, 356, 333, 590 和 242 个 氨基 
酸 。18S rRNA 的 基因 长 度 为 2 115 bp。 松 墨 天 牛 
的 这 9 个 内 参 基 因 和 其 他 昆虫 相应 基因 的 核 苷 酸 序 
列 一 致 性 非常 高 , 均 高 于 74% ,部 分 基因 的 一 致 性 
几乎 达到 100% ,说 明 这 些 内 参 基 因 在 进化 过 程 中 
非常 保守 。 








R1 松 墨 天 牛 内 参 基因 在 NCBI 上 的 Blastn 结果 


Table 1 Blastn results of reference genes of Monochamus alternatus 


















































基因 名 称 GenBank 登录 号 基因 长 度 (bp) 比 对 物种 比 对 基因 号 E ff 比 对 长 度 一 致 性 (% ) 
Gene name GenBank accession number Gene size Hit. species Hit. accession E-value Hit length Identity 
松 墨 天 牛 
TUB KU521361 1 579 EU073050 0 1 579 99.7 
Monochamus alternatus 
RPS27A KU521362 632 . AJ783900 2.7E-150 471 85.6 
Biphyllus lunatus 
RPLAO KU521363 754 : 杰 拟 谷 次 XM. 970465 2.7 E-169 666 80.2 
Tribolium castaneum 

AK KU521364 1 475 ADAE XM. 966707 0 1 430 76.9 

T. castaneum 
MEM 

GAPDH KU521365 1 335 TRE ! KC161214 0 978 76.9 
Bemisia tabaci 

EFIA KU521366 2 897 x BURMA NM. 001298945 0 1285 74.9 

Drosophila melanogaster 

RPS3 KU521367 832 杰 拟 谷 次 XM_008194244 0 759 79.8 

T. castaneum 
: 松 墨 天 牛 

Actin KU521368 1 600 EU100951 0 1 361 99.2 
M. alternatus 

18S rRNA KU521369 2 115 EXE JN619131 0 1 898 99.4 





Pogonocherus hispidulus 





2.2 引物 扩 增 效率 及 特异 性 

本 试验 利用 梯度 稀释 的 cDNA 作为 模板 进行 
RT-qPCR 扩 增 后 绘制 每 个 内 参 基因 的 标准 曲线 ,得 
到 每 个 基因 定量 分 析 引 物 的 相关 参数 ( 玉 ) ( 表 2)。 
从 表 2 中 可 以 看 到 9 个 内 参 基因 的 线性 相关 系数 变 
化 范围 为 0.994 ~1.000, 表 明 cDNA 初始 模板 量 与 其 
相应 的 Ct 值 线性 关系 明确 ,得 到 的 线性 相关 方程 可 
靠 。 HR E = (105779 -1) x100 计算 得 到 的 扩 增 
效率 为 102.3% ~ 109. 4% , 均 大 于 10096 ,表明 RT- 





qPCR 反应 的 有 效 性 。 模 板 与 Ct 值 间 的 良好 线性 关 
系 和 反应 的 有 效 性 确保 了 定量 结果 准确 可 靠 。 
2.3 不 同样 品 中 内 参 基 因 的 表达 水 平 差异 

本 研究 一 共 使 用 了 幼虫 . 晴 和 成 虫 的 12 个 组 
2H , RT-qPCR 结果 显示 如 图 1。9 个 内 参 基 因 在 研 
究 所 用 组 织 中 的 Ct 值 范 围 为 6.5 ~27.3, 表 明 9 个 
内 参 基 因 的 表达 量 均 较 高 ,其 中 ,18S rRNA 的 表达 
量 最 高 ,平均 Ct 值 为 10.0; EFIA 的 表达 量 最 低 , 平 
均 Ct 值 为 24.1; TUB 和 GAPDH 的 平均 Ct 值 分 别 
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R2 本 实验 所 用 引物 序列 、 扩 增 效率 及 线性 相关 系数 


Table 2 Sequence, amplification efficiency and correlation coefficient of the primers used in our study 




















ELA B 正 向 引物 (5" -3° 反 向 引物 (5' -3') UM E uti: MERECE 

Gene name Forward primer Reverse primer (bp ) : nplification Correlation 

Amplicon size efficiency coefficient 
18S rRNA AGAAACGGCTACCACATCCA GCGGGACACAGATTCGACTA 253 107.1 0.997 
Actin CGAACAGGAAATGGCTACCG CGTCGCACTTCATGATGGAG 191 112.3 0.997 
AK TCCCGCTGGTGAATTTGTTG AGCCTGCAAGAAACGATCAC 250 103.8 1.000 
EFIA GACAGCTCCATTTTCCACGTTG TTCGGGCGCTTTTCCAAATG 76 108.2 1.000 
GAPDH TCGAACGCTTCATGCACAAC CCATCACGCCACAATTTTCC 155 103.1 1.000 
RPLIO GACAAACGTTTCAGCGGAAC TGCTGTTGATCGCCCAAAAC 83 106.7 0.998 
RPS3 ATGTGGAGCTAAAGGTTGCG AAACGTTATCGGGCAAAGGC 242 108.9 0. 999 
TUB AATCACCAACGCCTGTTTCG TGGCAATAGCAGCGTTTACG 131 109.4 0.994 
RPS27A CGCCGGAAAACAATTGGAAG GCATTCACGTCTCAATCTGTGG 229 105.7 0. 999 


xo] 


254 


Ct 值 Ct value 


154 


104 
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18S rRNA Actin 




















AK EFIA GAPDH RPL10 RPS3 
基因 Genes 


到 1 9 个 候选 内 参 基 因 在 松 墨 天 牛 所 有 组 织 中 的 表达 比较 


TUB RPS27A 

















Fig. 1 Expression comparison of nine candidate reference genes in all the tested tissues of Monochamus alternatus 


2318.1 4120.3, 18S rRNA 在 不 同样 品 间 的 表达 量 
差异 最 大 ,Ct 值 的 最 大 差异 为 10.2; GAPDH f$ TUB 
在 不 同样 品 间 的 表达 量 差异 最 小 ,Ct 值 的 最 大 差异 
仅 为 3.4。 不 过 , 除 18S rRNA 和 Actin 外 ,四 分 位 差 
( interquartile range, IQR) 分 别 为 5.5 和 2.5 ,其 余 各 
内 参 基因 的 IQR 值 的 范围 为 1.3 ~1.7 ,十 分 接近 。 
2.4 geNorm 软件 分 析 

geNorm 通过 计算 基因 表达 稳定 值 M, 对 所 有 候 
选 内 参 基因 的 表达 稳定 度 进 行 排序 ,M 值 越 小 表明 
所 选 的 内 参 基因 越 稳定 。 由 图 2 可 知 ,对 于 松 墨 天 
牛 不 同 发 育 阶段 ,4K 和 GAPDH 是 最 稳定 的 内 参 基 

















参 基因 ,7UB 和 RPS 也 是 稳定 的 内 参 基因 。 内 参 
基因 18S rRNA FI Actin 在 不 同 发 育 阶 段 之 间 不同 
性 别 之 间 和 同一 性 别 不 同 组 织 之 间 的 稳定 值 均 
较 大 。 

在 基因 表达 分 析 时 ,往往 需要 使 用 2 或 2 个 以 
上 的 内 参 基因 对 目的 基因 进行 校正 ,从 而 减少 试验 
误差 ,得 到 更 加 准确 的 结果 。geNorm 软件 通过 分 析 
内 参 基 因 的 配对 差异 值 Vn/n + 1 变化 能 确定 在 某 
种 条 件 下 所 需 内 参 基因 的 最 适 数 目 。 本 研究 中 , 除 
了 V8/9 值 外 ,其 余 Vn/n +1 值 均 在 0.1 左右 ,小 于 
程序 推荐 值 0. 15( 图 3) ,说 明 在 对 不 同 发 育 阶段 之 




















因 ,7UB 也 是 稳定 的 内 参 基 因 。 对 于 不 同性 别 而 
言 ,TUB 和 GAPDH 是 稳定 的 内 参 基 因 。 对 于 同一 
性 别 的 不 同 组 织 而 言 ,4K 和 GAPDH 是 最 稳定 的 内 

















间 不 同性 别 之 间 和 同一 性 别 不 同 组 织 之 间 进 行 基 
因 表达 分 析 时 无 需 加 入 第 3 个 内 参 基因 进行 数据 
校正 ,最 合适 的 内 参 基 因数 目 是 2 个 。 
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图 2 GeNorm 软件 分 析 的 依次 删除 最 不 稳定 内 参 基因 后 所 有 基因 在 松 墨 天 牛 不 同 发 育 阶段 (A) 、 




















Hn B) 和 组 织 (C) 中 的 表达 稳定 值 M 


T T 1 
Actin RPS27A 18SrRNA AK EFIA RPL10 RPS3 GAPDH 


Fig. 2 Expression stability value M of all remaining reference genes after stepwise exclusion of the least stable reference genes 


右边 为 最 稳定 内 参 基因 , Z5 BU N fec E S HH 


in different developmental stages ( A) , genders (B) and tissues (C) of Monochamus alternatus 
因 。 虚 线 表示 表达 稳定 值 M 0.5, More stably expressed genes are positioned on the right 

















side, while less stably expressed on the left side. The dashed line indicates that the stability value M is 0. 5. 


Fig. 3 


2.5 NormFinder 软件 分 析 


成 对 变化 值 (Vn/n+1) 


Pairwise variation (Vn/n +1) between normalization factors after stepwise inclusion of stable reference 


Pairwise variation 





0.4 





四 发 育 阶段 Developmental stage O 性 别 Gender 8 组 织 Tissue 


0.3 


0.2 


0.1 





V2/3 V3/4 V4/5 V5/6 V6/7 V7/8 V8/9 
依次 增加 的 内 参 基 因数 量 


Number of reference genes 


图 3 从 最 稳定 内 参 基 因 开始 依次 增加 内 参 基 因 时 正 态 因子 间 的 成 对 变化 值 

















genes from the most stably expressed genes 
虚线 表示 成 对 变化 值 为 0. 15。The dashed line indicates that the pariwise variation is 0. 15. 


内 的 表达 稳定 值 来 筛选 内 参 基 因 , 以 筛选 最 佳 内 























NormFinder 通过 综合 计算 内 参 基 因 在 同一 组 。 参 基 因 。 稳 定 值 越 大 ,内 参 基因 越 不 稳定 ,稳定 值 
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越 小 ,内 参 基 因 越 稳定 。 本 研究 发 现 ,对 于 不 同 发 
育 阶段 和 不 同性 别 而 言 ,C4PD 的 稳定 值 最 小 ,对 
于 同一 性 别 的 不 同 组 织 而 言 , 7UB 的 稳定 值 最 
小 。18S rRNA 和 hctin 在 不 同 发 育 阶 段 . 不 同性 
别 和 同一 性 别 不 同 组 织 之 间 的 稳定 值 均 最 大 。 
奇怪 的 是 ,除了 不 同性 别 之 间 外 , NormFinder 推 
荐 的 用 于 不 同 发 育 阶 段 之 间 和 同一 性 别 不 同 组 
织 之 间 的 最 佳 引 物 组 合 均 不 含有 推荐 的 最 佳 内 
参 基因 ( 表 3)。 

表 3 通过 NormFinder 软件 分 析 获 得 的 候选 内 参 基因 在 

松 墨 天 牛 不 同 发 育 阶段 .性 别 和 组 织 中 的 表达 稳定 值 排序 

Table 3 The ranking order of the expression stability 








value of candidate reference genes in different 
developmental stages, genders and tissues of Monochamus 


alternatus calculated by NormFinder 





RIK FAE IÉ Expression stability value 











基因 








发 育 阶段 








Genes Developrmertil 性 别 ad 
ds Gender Tissue 

18S rRNA 0.221 0.320 0. 224 
Actin 0. 119 0.252 0. 069 
AK 0. 039 0. 161 0.041 
EFIA 0. 052 0. 163 0.037 
GAPDH 0.017 0. 103 0. 034 
RPL10 0. 065 0. 208 0. 038 
RPS3 0. 043 0. 140 0. 054 
TUB 0. 042 0.151 0. 024 
RPS27A 0. 089 0.250 0. 061 
最 佳 组 合 0. 025 0.070 0.016 


Best combination. (EFIA, TUB) (GAPDH, RPS3) (EFIA, RPL10) 








最 后 一 行 括号 中 为 NomFinder 软件 推荐 的 2 个 最 佳 组 合 内 参 基 因 。 











Two best reference genes were listed in the parentheses of the last row. 


2.6 BestKeeper 软件 分 析 

BestKeeper 软件 分 析 主 要 以 标准 变异 系数 
(SD) 变异 相 关系 数 (CV) .皮尔 逊 相关 系数 (CC) 
和 概率 值 (P) 来 判断 内 参 基因 的 表达 稳定 性 ,SD 值 
和 CV (ERR, CC 值 越 大 ,表明 该 基因 的 表达 越 稳 
定 。 从 表 4 中 数据 可 以 看 出 ,对 于 不 同 发 育 阶 段 而 
言 , 18S rRNA，Actin，EF1A4，RPL10，RPS3 和 
RPS274 由 于 其 标准 变异 系数 大 于 1, 并 不 适合 作为 
内 参 基 因 。 对 于 不 同性 别 而 言 ,Actin 和 EFIA 由 于 
其 标准 变异 系数 大 于 1,RPL10 由 于 其 概率 值 大 于 
0. 05 ,并 不 适合 作为 内 参 基 因 。 对 于 同一 性 别 的 不 
同 组 织 而 言 ,18S rRNA, Actin, AK 4 EFIA 由 于 其 
标准 变异 系数 大 于 1,18S rRNA, Actin, EFIA, 
RPL10, RPS3 和 RPS274 由 于 其 标准 变异 系数 大 于 
1 ,RPL10 和 RPS3 由 于 其 概率 值 大 于 0. 05 , 均 不 适 
合作 为 内 参 基因 。 
2.7 以 GAPDH 和 TUB 为 内 参 时 各 基因 的 表达 量 
比较 

分 别 利 用 GAPDH 和 TUB 为 内 参 ,计算 各 基因 
在 各 组 织 中 的 相对 表达 量 (图 4)。 以 GAPDH 为 内 
参 时 ,7TUB 在 各 组 织 中 的 相对 表达 量 差异 较 小 ,最 
大 差异 4.7 倍 。RPS3 在 各 组 织 中 的 相对 表达 量 差 
异 也 较 小 ,最 大 差异 5. 1 倍 。 其 余 各 基因 在 各 组 织 
中 的 相对 表达 量 差异 均 较 大 ,最 大 差异 倍数 均 大 于 
11 fj; LA TUB 为 内 参 时 ,G4PDH 在 各 组 织 中 的 相 
对 表达 量 差异 较 小 ,最 大 差异 4.8 倍 。 其 余 各 基因 
在 各 组 织 中 的 相对 表达 量 差异 均 较 大 ,最 大 差异 倍 
数 均 大 于 8.7 倍 。18S rRNA 和 Actin 在 各 组 织 中 的 
相对 表达 量 差异 很 大 ,最 大 差异 倍数 均 大 于 100 倍 。 
























































表 4 利用 BestKeeper 软件 评价 9 种 内 参 候选 基因 在 松 墨 天 牛 化 学 感受 组 织 中 的 稳定 性 


Table4 Stability of nine candidate reference genes in the olfactory tissues of Monochamus alternatus based 


on quantification cycle values ( Cq) evaluated using the BestKeeper software 

















188 rRNA. Actin AK EFIA GAPDH RPLIO RPS3 TUB ~ RPS27A 
SD [ + CP] 2.59 1.43 0.99 1.84 0.82 1.19 1.11 0. 63 1.45 
发 育 阶 段 CV [% CP] 24.89 8.11 4.42 — 7.60 | 4.06 — 5.58 | 5.96 3.48 7.56 
Developmental stage CC [r] 0. 63 -0.23 0.65 093 0.89 0.02 0.53 0.41 -0.29 
P 0. 00 0.20 | 0.00 | 06.00 0.00 0.93 0.00 0.02 0.09 
SD [ + CP] 0.46 1. 25 0.9 111 — 0.99 0.588 0.8 — 0.72 0.69 
性 别 CV [% CP] 5.87 7.11 4.00 473 4.84 2.2 4.50 3.94 3.46 
Gender CC [r] 0.78 0.55 070 | 0.88 0.89 0.32 0.69 0.82 -0.39 
P 0. 00 0.00 | 0.00 | 06.00 0.00 0.06 0.00 0.00 0.02 
SD [ + CP] 2.74 1.03 1.03 1.22 0.68 058 0.72 0.72 0.71 
组 织 CV [% CP] 28.37 6.10 449 5.10 — 3.34 2.62 3.83 3.89 3.53 
Tissue CC [r] 0.74 0.08 — 0.38 . 0.86 | 0.56 | 0.36 0.01 0.87  -0.39 
P 0. 00 0.66 0.03. 0.00 0.00 | 0.05 0.94 0.00 0.03 
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图 4 以 C4PDE(A) 和 TUB(B) 为 内 参 时 各 基因 在 松 墨 天 牛 不 同 组 织 中 的 相对 表达 量 


Fig. 4 Relative expression level of genes in different tissues of Monochamus alternatus relative to that of 


GAPDH (A) and TUB (B) as reference genes 

















基因 在 FAA 中 的 表达 量 为 对 照 , 定 为 1 ,在 其 余 组 织 中 的 表达 量 以 FAA 为 对 照 进行 计算 ;误差 线 表 示 标 准 差 。The expression level in female 


adult antenna was arbitrarily defined as 1, and those in other tissues were normalized to them. Error bars represent standard error ( SE). FAA ; WEHE 




















触角 Antenna of female adult; FAL; WE, Hid Leg of female adult; FAM ; MEE FEE ( F5) Labium (maxilla) of female adult; LA; 幼虫 触角 


Larval antenna; LM; 4%) E F JE ( F) Larval labium ( maxilla); MAA ; 雄 成 虫 触 





f| Antenna of male adult; MAH ; 雄 成 虫 头 部 Head of male adult; 














MAL; HEREJE Leg of male adult; MAM: fi, Hz TF JE ( F) Labium ( maxilla) of male adult; MAS; 雄 成 虫 腹部 末节 The last abdominal segment 
of male adult; MPA; 雄师 触角 Male pupal antenna; MPM.: Afi FJEE( F5.) Male pupal labium (maxilla). 


3 讨论 


目前 常用 的 RT-qPCR 内 参 基因 包括 GAPDH, 
actin, RPS27A , i7 3& i & i ( ubiquitin conjugating 
UBC ) 、18S rRNA,28S 核糖 体 RNA 
(28S rRNA) , EFIA, TUB, i YE $E K D] F 4E 
(translation. elongation factor gene, TEF) , RPLIO, 
AK ,RPS3 等 ( Nolan et au/.，2006 ) 。 研 究 发 现 , 内 人 参 
基因 在 不 同 物种 ,不 同 发 育 阶段 和 不 同 组 织 中 的 表 


enzyme gene, 





达 稳 定性 存在 差异 ,应 该 根据 研究 物种 和 研究 目的 
选择 合适 的 内 参 基 因 (de Jonge et al., 2007; Sun et 


al., 2009 ; Huis et al., 2010) 。 松 墨 天 牛 虽 然 是 林业 
上 的 重大 害虫 ,但 是 GenBank 中 仅 能 找到 TUB 
2008) 和 Actin 2 个 常用 内 参 基 因 的 序 
列 。 我 们 在 前 期 转录 组 测序 结果 中 ,新 筛选 到 
RPS27A, RPL10, AK, GAPDH, EFIA, RPS3 和 18S 
rRNA 7 个 常用 内 参 基因 ,前 6 个 编码 氨基 酸 的 基因 
均 包 含 完整 的 开放 阅读 框 。 松 墨 天 牛 的 9 个 常用 内 
参 基 因 和 其 他 昆虫 相应 基因 的 相似 性 非常 高 ,说 明 
这 些 基因 在 进化 过 程 中 非常 保守 。 

松 墨 天 牛 基因 表达 研究 中 , 使 用 的 内 参 基 因为 





(Song et al., 








B-actin (Lin et al., 2009; Xu et al., 2014; Luo et al., 
2015) fll 18S rRNA (Qian et al., 2015) 。 肌 动 蛋 白 
是 真 核 生物 中 广泛 存在 的 一 种 蛋白 质 ,是 微 丝 的 结 
构成 分 ,在 真 核 细 胞 进化 过 程 中 高 度 保守 ,被 广泛 用 
于 基因 表达 的 内 参 基 因 (Rodrigues ef al., 2014; 
Kumari et al., 2015) 。 核 糖 体 RNA 是 核糖 体 中 起 主 
要 作用 的 结构 成 分 ,存在 于 所 有 生物 中 ,被 认为 是 细 
胞 中 最 稳定 表达 的 基因 (Smit et al., 2007) ,在 研究 
人 类 组 织 细胞 时 , 常 被 作为 内 参 基 因 用 于 基因 表达 
人 研究 ( Goidin et al., 2001), KIR EM, actin. 和 
18S rRNA 在 各 化 学 感受 组 织 间 的 表达 差异 均 大 于 

其 他 候选 内 参 基 因 ，geNom，NormFinder 和 
BestKeeper 的 分 析 结 果 表 明 B-actin 和 188 rRNA 在 
所 用 化 学 感受 组 织 间 的 表达 非常 不 稳定 ,否定 了 B- 
actin 和 18S rRNA 作为 松 墨 天 牛 化 学 感受 组 织 基因 
表达 研究 的 内 参 基 因 。Actin 和 18S rRNA 不 适合 作 
为 基因 表达 研究 的 内 参 基因 在 其 他 物种 中 也 有 报 
道 , 例 如 actin 不 适合 作为 家 看 Bombyx mori, 甜菜 夜 
WS. exigua 和 LESSER mL. decemlineata [8] A &j 
阶段 的 基因 表达 内 参 基 因 ( Teng et al., 2012; Shi et 
al., 2013) ,18S rRNA 不 适合 作为 Bt 毒素 诱导 小 菜 
蛾 基因 表达 研究 的 内 参 基因 (Fu et al., 2012) 。 
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GAPDH 是 生物 体内 糖 酵 解 糖 异 生 及 光合 作用 
碳 固定 循环 过 程 中 的 关键 酶 ,在 进化 上 十 分 保守 。 








感受 组 织 中 的 内 参 基 因 进 行 筛选 , 发现 GAPDH 
松 墨 天 牛 化 学 感受 基因 表达 分 析 的 首选 内 参 基 





本 研究 发 现 ,使 用 的 9 个 候选 内 参 基因 中 ,GA4PDH 
在 不 同样 品 间 的 Ct 值 差 异 最 小 , geNorm 和 
NormFinder 软件 计算 的 基因 表达 稳定 值 最 低 ,同时 
完全 符合 BestKeeper 规定 的 作为 内 参 基因 的 要 求 ， 
说 明 C4PD 厅 在 松 墨 天 牛 不 同 发 育 阶段 不同 性别 
的 化 学 感受 组 织 中 的 表达 非常 稳定 ,是 松 墨 天 牛 化 
学 感受 组 织 表达 研究 中 最 合适 的 内 参 基 因 , 可 以 用 
于 不 同 发 育 阶段 .不 同性 别 和 同一 性 别 不 同 组 织 的 
化 学 感受 组 织 中 基因 表达 研究 。G4PDH 的 组 织 
达 稳 定性 也 在 B. mori 中 得 到 了 证 实 (Teng et al., 
2012) 。 虽 然 有 研究 发 现 GAPDH TE TENE SI T] 28 2H 
中 的 表达 不 稳定 (Jacob et al., 2013) ,但 也 有 研究 指 
出 GAPDH 在 不 同类 型 的 人 类 胎盘 中 的 表达 非常 稳 
^E ( Murthi et al., 2008) ,造成 这 种 差异 的 原因 可 能 
是 由 于 疹 椎 动物 基因 组 中 存在 大 量 的 GAPDH 假 基 
( pseudogene ) (Sun et al., 2012) 。 

TUB 与 细胞 形态 的 维持 、 细 胞 内 物质 的 运输 、 
细胞 运动 及 细胞 分 裂 增殖 等 生命 活动 密切 相关 。 本 
研究 发 现 ,TUB 基因 在 不 同样 品 间 的 表达 差异 也 很 
小 ,geNorm 计算 发 现 TUB 在 不 同 发 育 阶段 .不 同性 
别 和 同一 性 别 不 同 组 织 间 的 基因 表达 稳定 值 均 低 于 
0. 5, NormFinder 软件 计算 的 表达 稳定 值 也 较 低 , 完 
全 符合 BestKeeper 规定 的 作为 内 参 基因 的 要 求 ， 
此 TUB 也 是 松 墨 天 牛 化 学 感受 组 织 基因 表达 研究 
中 的 优秀 内 参 基 因 。7UB 作为 基因 表达 内 参 基因 
也 在 长 红 锥 晴 Rhodnius prolixus (Paim et al., 2012) 
和 棉花 粉 蛤 Phenacoccus solenopsis ( Chen and Lu, 
2014) 中 得 到 证 实 。 

使 用 单一 内 参 基因 进行 RT-qPCR 数据 校准 和 
标准 化 ,可 能 会 对 结果 的 精确 性 产生 影响 , 因此 在 
RT-qPCR 数据 标准 化 过 程 中 ,应 该 使 用 多 个 内 参 基 
因 进 行 校正 (Vandesompele et al., 2002) 。 本 人 研究 发 
现 , 当 增 加 内 参 基因 数量 时 ,geNorm 软件 计算 的 配 
对 差异 值 变化 小 于 程序 推荐 值 0.15 ,因此 在 研究 松 
墨 天 牛 化 学 感受 组 织 基 因 表 达 时 ,最 合适 的 内 参 基 
因数 目 是 2 Mo GAPDH 和 TUB 是 松 墨 天 牛 化 学 感 
受 组 织 中 表达 最 稳定 的 两 个 内 参 基因 ,其 表达 水 平 
相近 ,满足 内 参 基 因 组 合 时 的 各 项 要 求 ( Yuan et al., 
2012) ,是 松 墨 天 牛 化 学 感受 组 织 基 因 表 达 分 析 时 
的 推荐 内 参 基 因 组 合 。 

本 研究 从 松 墨 天 牛 中 新 鉴定 出 7 个 候选 内 参 基 
,并 利用 RT-qPCR 技术 和 3 种 软件 首次 对 不 同化 

























































































Dy pm 4k 


,内 参 基 因 G4PDH fü TUB 组 合 足以 对 松 墨 天 牛 
化 学 感受 基因 的 表达 进行 有 效 校正 。 由 于 内 参 基因 
的 稳定 性 与 物种 类 别 、 实 验 条 件 、 组 织 类 型 等 有 很 大 
相关 性 ( Sun et al., 2009; Huis et al., 2010) ,本 研 
究 基 于 松 墨 天 牛 化 学 感受 组 织 进行 的 内 参 基 因 筛 选 
结果 ,是 否 适用 于 其 他 昆虫 化 学 感受 组 织 以 及 松 墨 
天 牛 其 他 组 织 的 RT-qPCR ,还 有 待 于 进一步 研究 。 
在 进行 RT-qPCR 内 参 基因 选择 时 ,应 根据 特定 的 研 
究 目的 仔细 分 析 和 评价 内 参 基因 的 稳定 性 (de 
Jonge et al., 2007 ) 。 本 人 研究 结果 为 利用 RT-qPCR 
技术 准确 分 析 松 墨 天 牛 和 其 他 天 牛 基因 包括 化 学 感 
受 组 织 基因 相对 表达 量 的 内 参 基 因 选 择 提 供 了 重要 
的 参考 。 
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